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William A. Wood* Christopher .7. Riley 1 and F. McNeil Cheatwood* 
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Three equilibrium-air numerical solutions are presented for the Reentry-F flight-test 
vehicle at Mach 20, 80,000 ft. conditions, including turbulent flow predictions. The three 
solutions are from a thin-layer Navier-Stokes code, coupled thin-layer and parabolized 
Navier-Stokes codes, and an approximate viscous shock-layer code. Boundary-layer and 
shock-layer profiles are presented and compared between the three solutions, revealing 
close agreement between the three solution methods. Notable exceptions to the close 
agreement, with 7-10 percent discrepancies, occur in the density profiles at the boundary- 
layer edge, in the boundary-layer velocity profiles, and in the shock-layer profiles in 
regions influenced by the nose bluntness. 


Nomenclature 

c Sound speed, m/s 

M Mach number 

P Pressure, Pa, 

Re s Reynolds number based on length s 

T Temperature, K 

V Velocity, m/s 

j\ y. z Cartesian body axes, m 

a Angle of attack, deg 

// Wall-normal distance, m 

p Density, kg/m 3 

Subscripts: 

n Normal 

t Tangential 

tran Transition 

w Wall 

dc Freest ream 

Introduction 

R EENTRY-F 1 was a flight test conducted in 1968 
to study the turbulent reentry environment. The 
reentry vehicle itself consisted of an instrumented 5- 
deg half-angle cone, 3.66 m (12 ft.) in length. 

The 25 years since the original test have seen a 
tremendous advancement in the computational re- 
sources available for data analysis. The present study 
seeks to revisit the issues raised by the Reentry-F 
flight with regards to hypersonic laminar-to-turbulent 
transition by providing high-fidelity coinplete-flowfield 
solutions about the vehicle. These flowfield solutions 
can then provide a basis for the development of im- 
proved transition prediction and modeling methods. 

The application of such improved transition and 
turbulence prediction methods would be of signifi- 
cant benefit to the X-33 and X-34 reusable- launch- 
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vehicle programs. Current analysis techniques for hy- 
personic turbulent flow have appreciable aerodynamic 
and aerothermal uncertainties, and improved predic- 
tion methods would lower the costs associated with 
over-designing the vehicles to compensate for the pre- 
diction uncertainties. 

Codes 

Three algorithms, each representing a different level 
of physical modeling, were used to generate the flow- 
field solutions in the present study. 

The* Langley Aerothermodvnamic Upwind Relax- 
ation Algorithm (LAURA) 2 is a finite- volume, shock- 
capturing algorithm with second-order spatial accu- 
racy for the steady-state solution of viscous or inviscid 
hypersonic flows. The scheme employs a point implicit 
relaxation strategy with upwind-biased flux-difference 
splitting for perfect gas, equilibrium air, or nonequi- 
librium air calculations. In the present study LAURA 
solved the thin-layer Navier-Stokes (TLNS) equations 
for laminar and Far ve- averaged flows. 3 

Solving the TLNS equations is time and memory 
intensive, iterating globally on the domain in a time- 
marching fashion. However, this is a robust approach 
that is good for handling flowfields with subsonic re- 
gions. 

LAURA was also employed in conjunction with 
the Upwind Parabolized Navier-Stokes Solver (UPS). 3 
The procedure of Wood 5 " 8 was used whereby the flow- 
field about the spherically-blunted nose (3.43 mm ra- 
dius) was solved with LAURA while the 4 m conical 
body was solved with UPS. UPS is an upwind, finite- 
volume, parabolized Navier-Stokes (PNS) code with 
perfect gas, equilibrium air, and chemical nonequi- 
librium capability. It is second-order accurate in the 
crossflow plane and first order accurate in the march- 
ing direction, and contains a Baldwin-Lomax* alge- 
braic turbulence option. 

The PNS equations are cast as a space-marching 
scheme, solving the entire domain in a single stream- 
wise sweep. This is contrasted with a TLNS algorithm 
which typically iterates hundreds or thousands of times 
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on the domain. However, PNS algorithms are lim- 
ited to fully supersonic flows, with the exception of 
the boundary layer which must be treated specially 
in most PNS formulations. UPS uses the method 
of Vigneron 10 to stabilize the subsonic portion of the 
boundary layer. 

Since Reentry-F is modeled with a spherical nose 
in the present study, the PNS code cannot be utilized 
upstream of the sonic line. This is why LAURA is 
needed to solve the nose region before UPS can march 
the solution down the conical portion of the vehicle. 

The third code utilized was the Approximate Vis- 
cous Shock-Layer (AVSL) 11 algorithm. AVSL is an 
approximate viscous shock-layer technique whose gov- 
erning equations are identical to those of the standard 
viscous shock-layer technique except that Maslen’s 
pressure relation 12 is substituted for the normal mo- 
mentum equation. The streamwise momentum equa- 
tion remains unchanged. 

AVSL is a shock- fitting inverse method, iterating 
on the bow-shock shape to match the desired surface 
shape for that portion of the shock layer encompassed 
by the sonic line. Downstream of the sonic line the 
solution is space-marched in a manner similar to the 
PNS treatment, and hence requires a supersonic shock 
layer downstream of the nose. Being shock-fitting, 
AVSL is not appropriate to flowfields with internal 
shocks, whereas both UPS and LAURA, being shock- 
capturing, can model the internal shocks in addition 
to the bow shock. 

Cases 

Three cases are considered, corresponding to the 
three sets of governing equations solved by the 
LAURA, UPS. and AVSL computer codes. All of the 
solutions simulate the trajectory point at an altitude 
of 24.38 km (80,000 ft.). 

Case A is the three-dimensional combined LAURA- 
UPS solution. Tin 1 angle of attack is 0.14 deg. Both 
LAURA and UPS were run using the equilibrium air 
model of Srinivasan, Tannehill, and Weilmuenster. 13 * 
The UPS solution was begun at 3.55 mm and extended 
beyond the full vehicle length to 4 m. Transition was 
specified to commence at 2 m. in approximate agree- 
ment with the flight data, giving Re s = 36.9 x 10 6 
based on surface length from the nose. The wall 
temperature was held constant at 467.9 K. Table 1 
lists the freest ream conditions for Case A along with 
Cases B and C. The computational mesh contained 
101 points normal to the body surface and 15 points 
in the circumferential direction around half the vehicle. 
Only the port side of the vehicle was included in the 
computational domain since there was zero side-slip. 
The final streamwise step-size for the UPS portion of 
the solution was 5 111111 . 

Case B is an axisymnietric LAURA solution. In this 
case the equilibrium air model of Liu and Vinokur 14 


Table 1 

Reentry-F freestream conditions 


Case A 

Case B 

Case C 

M oc 

19.97 

19.97 

19.97 

a 

0.14 

0 

0 

Poo 

0.04479 

0.446 

0.04480 

Too 

221.3 

221 

221.3 

T w 

467.94 

variable 

467.94 

Poo 

2830.9 

2843.5 

2844.6 

Voo 

5936.4 

5966 

5954.2 


297.3 

298.7 

298.2 

%end 

4.036 

3.683 

4.036 

•E tran 

2.012 

2.0 

2.012 


Table 2 T u , variation for Case B 

X T u , 

0 

1453 

0.053 

1194 

0.078 

1065 

0.130 

806 

0.257 

667 

0.792 

513 

1.328 

462 

1.864 

437 

2.527 

427 

3.086 

486 

3.673 

542 


was used along with the Cebici-Smith 15 algebraic tur- 
bulence model. As in Case A, transition was specified 
to begin at 2 m. The wall temperature varied along 
the length of the vehicle according to the measured 
flight data with a maximum of 1450 K at the nose and 
a minimum of 427 K located just past the onset of 
transition. A tabulation of the imposed wall temper- 
ature versus axial distance for Case B is contained in 
Table 2. The axisymmetric grid contained 201 points 
in the streamwise direction and 97 points normal to 
the body. 

Case C is an axisymmetric AVSL solution. The equi- 
librium air model used is due to Hansen. 16 The Cebici- 
Smith turbulence model is employed with AVSL and 
the transition region is modeled after Dhawan and 
Narashima. 1 ' This is the same transition method as is 
used in both UPS and LAURA. The wall temperature 
was held constant at the same value as for Case A. 
The grid for this case contained 101 points normal to 
the body. 

Results 

Results will be presented in the form of profile plots 
taken from the windside and leeside centerlines of 
Case A, the three-dimensional case, and the corre- 
sponding points from Cases B and C, the axisymmetric 
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Fig. 1 Shock-layer density profiles, x = 0.5 m. 


UPS - windside 

— — - UPS - leeside 



Fig. 2 Shock-layer density profiles, x = 1.0 m. 

cases. Both shock- layer and boundary- layer profiles 
are displayed versus the body-normal distance from 
the surface. The profiles presented are for density, 
f>!() oc, surface- tangential velocity component, \)/c oc , 
surface-normal velocity component, 1 ,,/cx, and pres- 
sure, PjP 0 c ,. Profiles are extracted at five axial loca- 
tions, x — 0.5, 1.0, 1.5, 2.0, and 4.0 m. 

Density profiles 

Shock-layer density profiles are plotted in Figs. 1 5 
at the five axial data stations. Figure 5 plots only the 
results of Cases A and C at x = 4 in, because Case B, 
the LAURA solution, had been terminated prior to 
this station. The profiles from Cases A and B extend 
beyond the shock layer into the freestream because 
both UPS and LAURA are shock-capturing codes. 
The results of Case C stop at the shock because AVSL 
is a shock-fitting code, obtaining post-shock conditions 
analytically from the Rankine-Hugoiiiot relations. 

Figure 1 shows the bow shock is better resolved at 
the most upstream location, x = 0.5 in, by the LAUR A 
solution than the UPS result. This is attributed to a 


UPS - windside 
UPS - leeside 



Fig. 3 Shock-layer density profiles, x = 1.5 m. 


UPS - windside 


— — - UPS - leeside 



Fig. 4 Shock-layer density profiles, x = 2.0 m. 


UPS - windside 


— — - UPS - leeside 



Fig. 5 Shock-layer density profiles, x = 4.0 m. 
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Fig. 6 Boundary-layer density profiles, x = 0.5 m. 

tighter grid-spacing normal to the shock in the Case B 
versus Case A grid. Figures 2 4 show comparable 
shock-resolution between LAURA and UPS, where the 
grids are more similar than at x = 0.5 m. 

In Fig. 1 the shock-layer is still strongly influenced 
by the nose bluntness, as the density profile outside 
the boundary layer is continuously curved. Progress- 
ing downstream in Figures 2 5 the shock-layer assumes 
a more conical-flow profile, becoming nearly linear out- 
side the boundary layer. However, a “bubble” in the 
density profile does persist between the boundary layer 
and the conical-flow profile downstream to the last, 
laminar station, Fig. 4. This “bubble’ 1 persisting since 
the nose is swallowed by the turbulent boundary layer, 
Fig. 5. 

Comparing the performance of the three codes, 
LAURA predicts higher densities bv 7 percent than 
either UPS or AVSL, which agree closely with each 
other, at x = 0.5 m in Fig. 1. In all figures the wind- 
side and leeside results of Case A are very similar, 
which is not unexpected since the angle of attack is 
only 0.14 deg. The UPS and LAURA profiles over- 
plot in Figs. 2 4. The AVSL solution matches the 
UPS and LAURA profiles except that the “bubble” 
above the boundary layer is lost by x = 1 in, seen 
in Figs. 2 4. The turbulent profile* in Fig. 5 shows a 
3-percent difference in density between the UPS and 
AVSL solutions. 

Boundary-layer density profiles are charted in 
Figs. G 10, following the same sequence as Figs. 15. 
At the furthest upstream station, x = 0.5 m. Fig. G, 
excellent, agreement through the boundary layer is seen 
between the Case A UPS solution and the Case B 
LAURA result. At the boundary-layer edge the UPS 
and LAURA solutions separate, establishing the dif- 
ferent shock-layer profiles seen in Fig. 1. The Case C 
AVSL solution has excellent agreement throughout 
two-thirds of the boundary layer, but does not recover 
the density to the shock-layer value as soon as UPS 
and LAURA predict. The AVSL densities over-state 



Fig. 7 Boundary-layer density profiles, x = 1.0 m. 



Fig. 8 Boundary-layer density profiles, x = 1.5 m. 



Fig. 9 Boundary- layer density profiles, x = 2.0 m. 
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Fig. 10 Boundary-layer density profiles, 

x = 4.0 m. 



Fig. 12 Shock-layer tangential- velocity profile, 
x — 1.0 m. 



Fig. 11 Shock-layer tangential- velocity profile, 
x — 0.5 in. 

the influence of the boundary layer by 15 percent, a 
trend repeated in Figs. 7-9. 

Density profiles at the remaining laminar stations, 
Figs. 7-9, show excellent agreement between the three 
cases out to the boundary-layer edge, where the 
LAURA edge profile is sandwiched bv the UPS edge 
profile below and the AVSL profile above. The three 
cases converge to similar shock-layer values outside the 
boundary-layer edge region. 

At the turbulent station, Fig. 10, the AVSL solution 
predicts up to 10 percent higher densities than UPS 
through the boundary layer and up to the boundary- 
layer edge, when* the differences are larger. 

Tangential- velocity profiles 

Shock-layer profiles for the velocity component tan- 
gential to the cone surface, in the streamwise direction, 
are plotted in Figs. 1 1 15. These figures show excellent 
agreement between all three solution methods through 
the shock layer. This agreement contrasts somewhat 
with the differences noted in the shock-layer density 


0.04 r 


0.03 h 
rj, m 


UPS - windside 

— — - UPS - leeside 

LAURA 

AVSL 


0.01 h 


I I — T 


10 

V./c 


20 


Fig. 13 Shock-layer tangential- velocity profile, 
x — 1.5 m. 



Fig. 14 Shock-layer tangential- velocity profile, 
x = 2.0 m. 


5 of 1 1 


NASA Langley Research Center 




Fig. 15 Shock-layer tangential-velocity profile, 
x = 4.0 m. 



Fig. 16 Boundary- layer tangential- velocity profile, 
x = 0.5 m. 

profiles. The densities are more sensitive than the ve- 
locities to the differing equilibrium chemistry models 
in the three ('odes and to predicted temperature vari- 
ations throughout the shock layers. 

Boundary-layer profiles of the tangential velocity are 
presented in Figs. 16 20, revealing a noticeable differ- 
ence between the UPS solution and the LAURA and 
AYSL solutions. At each of the four laminar stations. 
Figs. 16 19, UPS predicts a thinner boundary layer, by 
10 percent at x = 2.0 m, with a fuller profile through- 
out the boundary layer. This contrasts somewhat with 
the density boundary-layer profiles in Figs. 6 9 where 
the densities for all three cases were* similar through 
the majority of the boundary layer, differing only in 
the edge location. 

Examining more closely as a typical location the 
x = 1.5 m station. Fig. 18, UPS predicts a 7 per- 
cent increase in tangential velocity through the outer 
two- thirds of the boundary layer, relative to either the 
Case B or Case C results. The three solution sets 
converge at the boundary- layer edge and at the wall, 

6 



Fig. 17 Boundary-layer tangential- velocity profile, 
x = 1.0 m. 



Fig. 18 Boundary-layer tangential- velocity profile, 
x = 1.5 m. 



Fig. 19 Boundary-layer tangential- velocity profile, 
x — 2.0 m. 
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Fig. 20 Boundary-layer tangential- velocity profile, 
x = 4.0 m. 



Fig. 21 Shock-layer normal- velocity profile, 
x — 0.5 m. 

where the no-slip condition constrains all of the solu- 
tions. Differences in the implementation of the equi- 
librium thermal-property models in the three codes 
could be contributing to the different velocity gradi- 
ents through the boundary layer, with a decrease in 
viscosity corresponding to an increase in tangential- 
velocity gradient. 

The turbulent tangential- velocity profiles, Fig. 20, 
differ between UPS and AVSL. The AVSL profiles 
is much fuller than the UPS profile, though the 
boundary-layer edge predictions are similar. 

Normal- velocity profiles 

Shock-layer profiles of the velocity component nor- 
mal to the vehicle surface are presented in Figs. 21-25. 
The shock-layer values of V„ are negative, indicating 
flow toward the vehicle 1 surface*. A linear variation in 
U„ is seen through the shock layer, consistent with 
conical flow. 

Figure 21, the furthest upstream location, shows a 
slight variation in the AVSL solution from the LAURA 



Fig. 22 Shock-layer normal- velocity profile, 
x = 1.0 m. 



Fig. 23 Shock-layer normal- velocity profile, 
x = 1.5 m. 



Fig. 24 Shock-layer normal- velocity profile, 
x — 2.0 m. 
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Fig. 25 Shock-layer normal- velocity profile, 
x = 4.0 ni. 



V p /C. 

Fig. 26 Boundary- layer normal-velocity profile, 
x = 0,5 m. 

results near the shock. Excellent agreement between 
LAURA and AVSL is seen at three remaining laminar 
stations. Figs. 22 24. 

A very slight , yet discernible, difference can be seen 
between the UPS and LAURA profiles in Figs. 21-24. 
Both Case A and B results show the same slope of 
1 „ vs //. but UPS predicts slightly larger magnitudes 
of the normal-velocity component at a given height 
above the surface. 

At the t urbulent station. Fig. 25, the normal velocity 
is positive near the surface, due to the rapidly thick- 
ening displacement thickness in the turbulent bound- 
ary layer. The AVSL and UPS solutions agree well 
throughout the shock layer in Fig. 25. 

Boundary-layer profiles for the normal- velocity com- 
ponent an 1 shown in Figs. 26-30. Normal velocities 
are seen to be nearly zero in the laminar boundary 
layers between x = 1-2 m, Figs. 27-29. UPS pre- 
dicts slight normal velocities in this region, whereas 
the LAURA and AVSL solutions remain at zero. This 
small difference in the UPS profile could be related to 



V n /C . 


Fig. 27 Boundary-layer normal- velocity profile, 
x — 1.0 m. 



V c 


Fig. 28 Boundary- layer normal-velocity profile, 
r = 1.5 m. 


0.004 - 



V n /c. 


Fig. 29 Boundary-layer normal- velocity profile, 
x = 2.0 m. 
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Fig. 31 Shock-layer pressure profile, x = 0.5 m. 



the imposition of the Yigneron condition on the sub- 
sonic portion of the UPS boundary layer. 

Further upstream at r = 0.5 in. Fig. 2G, the normal 
velocities in the boundary layer deviate linearly from 
zero, with the UPS magnitudes still larger than either 
LAURA or AVSL, which agree with eacli other. These 
observed non-zero boundary-layer normal-velocity val- 
ues are indicative of nose-bluntness influence at this 
station. Still, the magnitudes of the normal- velocity 
component in the boundary layer are two orders of 
magnitude smaller than the tangential-velocity com- 
ponents in Fig. 1G. 

The turbulent normal-velocity profile. Fig. 30, shows 
an inflection in the profile* with positive, out-flowing 
values in the outer boundary layer. The overall mag- 
nitudes, though, remain very near zero. 

Pressure profiles 

Shock- layer pressure profiles for Cases B and C are 
presented in Figs. 31 35. The Case A UPS pressures 
were unavailable for comparison. Excellent agreement 
through the shock layer between LAURA and AVSL 


Fig. 33 Shock-layer pressure profile, x = 1.5 m. 



P/P_ 

Fig. 34 Shock-layer pressure profile, x — 2.0 m. 
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Fig. 35 Shock-layer pressure profile, x = 4.0 m. 

is seen at the three downstream laminar stations, 
Figs. 32 34. The surface pressure remains constant 
along the body in this region at P/Poc = 5.75. 

Closer to the nose, Fig. 31. where nose-blunt ness ef- 
fects are still influential, there are differences between 
the AYSL and LAURA profiles. The AVSL surface 
pressure comes in 5 percent lower than the Case B 
LAURA solution. 

The turbulent pressure profile. Fig. 35, retains the 
same shape as the laminar profiles, but with a 9 per- 
cent higher surface pressure. The thicker turbulent 
boundary layer has pushed the shock further from the 
body, increasing the shock strength and compression 
ratio. 

Boundary-layer profiles of the pressure are not pre- 
sented because the pressures remained constant at 
their edge values for all stations investigated. 

Heating 

It is noted that surface heat-transfer rates from the 
Case A solution have been previously presented in 
Ref. 5. In that report the Case A heating rates are 
seen to compare very well over the laminar range with 
data derived from the flight test. Heating agreement 
between flight and computation is reasonable in the 
transitional domain and good in the turbulent region 
toward the aft of the vehicle. 

Summary of Results 

Three levels of physical modeling were applied to 
the solution of the Reentry-F flight-vehicle fiowfield. 
Conditions correspond to a flight altitude of 24.4 km 
(80,000 ft.) at Mach 20. Calculations were performed 
using equilibrium-air chemistry with algebraic turbu- 
lence models. Results are presented in the form of 
shock-layer and boundary-layer profiles of density, ve- 
locity components, and pressure. 

Nose-bluntness effects were observed at the 
x — 0.5 m station. 145 nose radii downstream of the 
vehicle tip, indicated by non-linear shock-layer pro- 


files about the conical vehicle. Densities predicted by 
LAURA were 7 percent larger through the shock layer 
than given by UPS or AVSL at this station. Further 
downstream, at x = 1-2 m, very close agreement is 
seen between the solution sets, particularly with the 
LAURA and AVSL solutions. Notable exceptions are 
for the boundary-layer-edge density profiles and the 
boundary-layer velocities, where the UPS values are 
7 percent greater than from LAURA or AVSL. 

Pressure profiles show a 5 percent discrepancy be- 
tween LAURA and AVSL at x = 0.5 m, but excellent 
agreement further downstream. Post-shock pressures 
are seen to rise 1 by 9 percent after the onset of transi- 
tion. 

The dataset presented herein is intended to pro- 
vide high-fidelity, state-of-the-art numerical solutions 
about the Reentry-F flight-test vehicle, which can be 
utilized in conjunction with the flight -test data to de- 
velop and improve transition-prediction models. The 
development and application of such models is consid- 
ered an important need for the efficient design of the 
X-33 and X-34 reusable launch vehicles. 
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